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We analyze within a semiclassical approximation the form factor for the fluctuations of quantum 
matrix elements around their classical average. We find two contributions: one is proportional to 
the form factor for the density of states, with an amplitude determined by the squared average of 
the matrix elements. The other is constant and related to the fluctuations of finite time classical 
trajectory segments around the phase space average. The results are illustrated for an observable 
in the quadratic Zeemann effect. 
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Much of the relation between classical chaos and signa- 
tures in quantum mechanics has been clarified by a com- 
bined study of semiclassical and random matrix theories 
[1,2]. Semiclassical analysis relates the mean of eigen- 
values and matrix elements to a classical phase space 
volume part, the Thomas- Fermi or Weyl term and the 
fluctuations around the mean to periodic orbit contribu- 
tions through the Gutzwiller trace formula and [3]. Ran- 
dom matrix theory focusses on the universal properties of 
the fluctuations and has predictions for Gaussian ensem- 
bles with orthogonal, unitary and symplectic symmetries, 
which are in good agreement with numerical observations 
on chaotic systems [2,4]. Berry [5] has connected both 
approaches in his analysis of the spectral form factor, 
the Fourier transform of the correlation function for the 
density of states. For short times, the form factor is dom- 
inated by isolated periodic orbits [6,7]. For intermediate 
and long times qualitative aspects of the form factor can 
be explained by semiclassical and quantum sum rules. A 
numerical demonstration of this relation will be included 
below. 

Our main purpose here is to analyze the form factor 
for matrix element fluctuations. Previous semiclassical 
analyses have focussed on the mean behaviour [8-10] and 
correlations in billiard wave functions [11], periodic orbit 
effects in wave functions [12,13] and periodic orbit corre- 
lations in fluctuation properties [14,15]. Random matrix 
theory predicts Gaussian fluctuations of individual ma- 
trix elements and thus Porter-Thomas distributions for 
transition strengths [2,4]. One question to be answered 
here is that on the semiclassical prediction for the form 
factor of a matrix element weighted density of states. 

Closely related is the problem of the size of quan- 
tum fluctuations around the classical Thomas-Fermi 
terms. Calculations based on this classical expression 
have proven rather succesful in molecular physics, but 
their accuracy must be taken as accidental since there 
have been no a priori estimates for the error. The present 
analysis will yield such an estimate and will quantify the 
expectation that the classical behaviour is recovered as 



?z — ^ . We will focus here on strongly chaotic systems, 
leaving the case of mixed phase space and further inves- 
tigations to a more detailed report [16]. 

We consider a spinless particle in a bounded domain, 
with quantum eigenvalues E n ordered in an increasing se- 
quence and eigenstates \n) , and some observable A, which 
will be assumed to be sufficiently smooth so that the 
classical limit A(p,q) exists and that the periodic orbit 
expressions of [15] can be applied. We will assume that 
the classical system is fully ergodic without any islands 
of stability. A fixed energy around a reference energy Eo 
will contain an increasing number of states as Eo —> oo, 
while the classical mechanics changes rather little for suff- 
ficiently small intervals. We then consider the density of 
states, weighted by the matrix elements and pro- 

jected onto an interval of width AE around energy Eo 
by a window function, 

PA(E) = J2( n \Mn)6(E- E n )w(E n ). (1) 

n 

For the theoretical analysis we use a Gaussian window 

w(E) = -=1 e -(E-E r/2AE- (2) 

V ' ^AE V ' 

the choice of window functions of course does not affect 
the final result but some of the prefactors in intermediate 
expressions. Our analysis concerns the auto-correlation 
function of the fluctuations of the density of states around 
its average ~p A and its Fourier transform, the form factor. 
The auto-correlation function, defined by 

C A (e) = JdE PA (E + e/2) PA (E -e/2)-p 2 A , (3) 

splits into a diagonal and an off-diagonal part, 

C A (c) = J2\{n\A\n)\ 2 w(E n f6(c) 

n 

+ {n\A\n){m\A\m)w{E n )w{E m )6{e - (E n - E m )) . (4) 
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A convenient normalization is to divide by the weighted 
number of states under the window, 



(5) 



where ~p is the mean density of states; the remaining fac- 
tors are specific for the Gaussian window. With this nor- 
malization, the strength of the delta function at the origin 
in the normalized function 



C A (e) = C A (e)/C N 



(6) 



becomes 1 for the density of states and the average of 
the square of the matrix elements (n|j4|n) in the general 
case. 

The form factor is the Fourier transform of the auto- 
correlation function, 



K(t) = / deC(e)exp(id/h) 



(7) 



Since the (^-function in C(e) is mapped into the large 
t behaviour of K{i), the density of states and the size 
of matrix elements can be read off from the large t be- 
haviour of the form factor. 

The semiclassical analysis starts from the Gutzwiller 
trace formula [3,17,15], 

PaME) = Pa + ^-J2 VpApSW-W* , (8) 

p 

where the first term is the smooth phase space average of 
the observable A and where the sum extends over all pe- 
riodic orbits and multiple traversals. This form assumes 
that the window AE around Eo is sufficiently small so 
that the actions can be expanded linearly in energy, giv- 
ing rise to the exponential in (8). The term 



A v 



dtA(p(t),q(t)) 



(9) 



denotes the integral of the observable along the periodic 
orbit p. The remaining term w p depends on the mon- 
odromy matrix, phases due to the action at energy Eo 
and the Maslov indices. The semiclassical approxima- 
tion for the form factor of the matrix element weighted 
density of states then is the absolute value squared of the 
Fourier transform of the product of (8) with the window 
function, divided by the normalization factor C'n- Upon 
division by the normalization C'n, and after cancellation 
of various prefactors related to the window function one 
finds within the diagonal approximation [5] 
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\w p \ 2 A 2 p 8 T (t-T p ). 



(10) 



where Th = 2irph is the Heisenberg time. 



The widths of the (^-function is related to the widths of 
the initial window, r = Ti/ AE. For short times t, the 6- 
functions do not overlap, and one can identify individual 
periodic orbits [6,7]. For longer times, the (^-functions 
overlap, forming a quasi-continuum. Since the ampli- 
tude of the orbits contributing to this quasi-continuum is 
given by purely classical quantities, one can evaluated its 
amplitude in hyperbolic systems under the assumption 
that the periodic orbits are densely distributed in phase 
space [18-21]. In particular, for the density of states with 
A = 1, the integrals A p become the periods T p and one 
can use the classical sum rule [18-21] 



J2\w p \ 2 T^r(t-T p ) 



(T P 2 ) 



t 



to obtain 



Ki, 



gt/T H 



(11) 



(12) 



where the factor g = 2 for orthogonal systems and g = 1 
for unitary systems is related to the presence or absence 
of unitary symmetries [5]. 

If the observable is not the identity, the result depends 
on the fluctuation properties of the classical integrals 
A p . Since long trajectories will ergodically sample phase 
space, one can expect that on average A p ~ AT p , where 
A denotes the classical ergodic average of the observable 
A. Assuming that correlations along trajectories decay 
sufficiently rapidly, as in hyperbolic systems, the distri- 
bution around this average will be Gaussian, 

P T (A)dA = Prob {A p £ [A, A + dA] and T p near T} 



\/2iraT 



exp -(A - AT) 2 /aT , 



(13) 



where the variance increases linearly with the period T 
of the orbits. Thus the average of the squares of A p near 
a period T becomes 



(A^) ~ A 2 T 2 + aT . 



(14) 



Combining this with the classical periodic orbit sum rule 
(11), one arrives at the diagonal approximation for the 
form factor 

K A ,.c(t) ~ g~A 2 t/T H + ga/T H (t < T H ) . (15) 

Thus the form factor for matrix elements contains in ad- 
dition to the linearly increasing part a constant off-set. 
This off-set is determined by the fluctuation properties of 
averages of the observable computed from classical tra- 
jectory segment. 

For longer times, one can appeal to the behaviour of 
the quantum correlation function [5] and argue that K(t) 
has to settle at a constant, related to the strength of the 
delta function in C(e). In a simple approximation one 
can assume that the semiclassical form is correct up to 
a time t c ~ Th / g and that K(t) keeps the value thus 
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attained also for larger t. For the density of states form 
factor, this assumption yields the correct expression for 
the Gaussian unitary ensemble. Since uniformly hyper- 
bolic systems seem to fall into the universality class of 
Gaussian random matrix models [2], we will below com- 
pare our results to the functional form of the random 
matrix theory form factor. 

The asymptotic amplitude relates matrix elements and 
classical phase space averages, 



(n\A\n) 2 ~ g A + ga/T H 



(16) 



The first term is to be expected from the classical phase 
space average. The second term is a correction present 
for finite densities of states (finite Th) only, and can be 
calculated from the fluctuations of classical trajectory 
segments. If the average of the matrix elements is sub- 
tracted, only this term survives and it then measures the 
variance of the quantum matrix elements. It varies like 
1/Th and vanishes as k — ► 0. In non-hyperbolic systems, 
the dependence on Th may be different, as explained in 

We test these predictions for eigenstates for hydrogen 
in a magnetic field (for reviews, see [22,23]). The peri- 
odic orbit contributions to the density of states have been 
studied in [15]. Due to scaling, the classical dynamics de- 
pends on a combination of energy E and scaled magnetic 
field 7 = B/B with B = 2.35 • 10 5 T only. At scaled en- 
ergy e = Ey~ 2 / 3 = —0.1, a large fraction of phase space 
is chaotic and the small elliptic islands that will still be 
present are not felt on the time scale relevant for the 
present analysis. In this sense, the system is hyperbolic. 
A numerical solution of the Schrodinger equation for 
quantizing magnetic field strengths z n = j n 1/ ' 3 in semi- 
parabolic coordinates produced 4996 converged eigen- 
states and matrix elements for the observable A = l/2r, 
the potential energy. For the numerical determination of 
the form factor, these eigenstates were divided into five 
lots of 1000 states each, modulated by a Hamming win- 
dow [24]. The Fourier transforms were calculated in the 
density of states as a functions of the scaling variable z, 
the corresponding Heisenberg time being Th = 2tt/Az, 
where Az is the mean spacing of the eigenvalues. 

Fig. 1 shows the form factor of the density of states 
and compares it to the form factor for the Gaussian or- 
thogonal ensemble from random matrix theory, 



K, 



GOE 



2t-tln(l + 2t) (t < 1) 

■tln((2t+l)/(2t-l)) (t > 1) 



(17) 



To eliminate the large fluctuations evident in the in- 
set, the quotient of the calculated Fourier transform 
and the GOE form factor was averaged over an inter- 
val [max{0, t — 0.6}, min{t max , t + 0.6}], weighted with a 
Parzen window [24]. The fitted curve thus is 



Kcalcjt) 
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The remaining fluctuations indicate the size of the devi- 
ations to be expected. 

The form factor for the matrix element weighted den- 
sity of states is shown in Fig. 2. Since the average of 
the matrix elements is non-zero, one can still see the 
linear increase. The constant term to be discussed be- 
low is smaller than the fluctuations, so that the overall 
form is again compatible with that of a GOE form factor. 
The asymptotic amplitude (0.068) compares favourably 
with the square of the average matrix elements (0.26 2 = 
0.0676). 

To highlight the second contribution to (16) from the 
fluctuations of the classical trajectory segments, we show 
in Fig. 3 the form factor for the fluctuations of the marix 
elements minus their average. Within the fluctuations, 
this function is constant, as expected. As the quantum 
numbers and Th increase, the constant decreases. For 
the Heisenberg times relevant to Fig. 3 there are still 
some correlations and the rate of decay is not 1/Th, but 
the size of the fluctuations is still given by the variance 
of classical trajectory segments [16]. 

In summary, a semiclassical analysis of the form factor 
of the matrix element weighted density of states has re- 
vealed a dominant term due to the phase space average 
of the classical observable and a second contribution due 
to the fluctuations of classical trajectory segments. The 
shape of the first term is found to be qualitatively simi- 
lar to form factor for the density of states (semiclassical 
theory does at present not allow one to go further). The 
second part opens the possibility of estimate from clas- 
sical trajectories not just the average of matrix elements 
but also the size of the quantum fluctuations around this 
average, with important consequences for the classical 
trajectory calculations in atomic and molecular physics 
[25]. Further discussions and examples, as well as an 
extension to non-hyperbolic systems, will be discussed 
elsewhere [16]. 
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FIG. 1. Spectral form factor K\(t) for the eigenvalues of 
hydrogen in a magnetic field at scaled energy e = —0.1. The 
different curves correspond to the eigenvalues 1 — 1000 (dot- 
ted), 1001 - 2000 (dash-dotted), 2001 - 3000 (long dashed), 
3001 -4000 (short dashes) and 3997-4996 (continuous). For 
comparison, a fit to the GOE form factor as determined from 
(17) is also shown (thick line). The inset shows the original, 
fluctuating data for the eigenvalues n = 3997 . . . 4996. 

FIG. 2. Form factor Ka(1) for the eigenvalues of hydro- 
gen in a magnetic field weighted with the matrix elements of 
A = l/2r. The GOE form factor shown for comparision was 
adjusted by the asymptotic amplitude. Otherwise same as 
Fig. 1. 

FIG. 3. Form factor Ka(1) for the eigenvalues of hydro- 
gen in a magnetic field weighted with the matrix elements of 
A = l/2r — A. Since now the average of the matrix elements 
vanishes, there is no linear increase, the form factor is con- 
stant. The original data are shown in the inset, the averages 
were computed with a Parzen window. 
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